Filler-Induced Composition Waves in Phase-Separating Polymer Blends 
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The influence of immobile filler particles (spheres, fibers, platelets) on polymer blend phase separa- 
tion is investigated computationally using a generalization of the Cahn-Hilliard-Cook (CHC) model. 
Simulation shows that the selective affinity of one of the polymers for the filler surface leads to the 
development of concentration waves about the filler particles at an early stage of phase separation 
in near critical composition blends. These "target" patterns are overtaken in late stage phase sepa- 
ration by a growing "background" spinodal pattern characteristic of blends without filler particles. 
The linearized CHC model is used to estimate the number of composition oscillations emanating 
from isolated filler particles. In far-off-critical composition blends, an "encapsulation layer" grows 
at the surface of the filler rather than a target pattern. The results of these simulations compare 
favorably with experiments on filled phase separating blend films. 

PACS: 64.75.+g, 68.55.Jk, 47.54.+r, 61.41.+e 



I. INTRODUCTION 

The bulk properties of misciblc fluid mixtures are char- 
acteristically insensitive to their microscopic fluid struc- 
ture near the critical point for phase separation, where 
the properties are governed by large scale fluctuations in 
the local fluid composition. Composition fluctuations oc- 
cur similarly for most near-critical fluid mixtures, so that 
the properties of these fluids are subject to a "univer- 
sal" description. This accounts for the success of simple 
mathematical models of critical phenomena (e.g., Ising 
model, </> 4 -field theory) that contain the minimal physics 
of these fluctuation processes. Although mixtures near 
their critical point are susceptible to external perturba- 
tions, the influence of microscopic heterogeneities tend to 
become "washed out" in the large scale fluid properties, 
apart from changes in critical parameters describing the 
average properties of the fluid (e.g., critical temperature 
and composition, apparent critical exponents, etc.). This 
situation changes, however, when the fluid mixture enters 
the two-phase region. The fluid is then far from equilib- 
rium and perturbations can grow to have a large-scale 
influence on the phase separation morphology. Pertur- 
bations in these unstable fluids can be amplified rather 
than washed out at larger length scales. Inevitably, the 
theoretical description of this kind of self-organization 
process is complicated by various non-universal phenom- 
ena associated with the details of the particular model 
or experiment under investigation |l|. The beneficial as- 
pect of this sensitivity of phase separation and other pat- 
tern formation processes to perturbations is that it offers 
substantial opportunities to control the morphology of 
the evolving patterns and leads to a great multiplicity of 
microstructures. 

Many previous studies have considered the application 
of external influences (flow g and gravitational || fields, 



concentration J I| and temperature || gradients, chemi- 
cal reactions [| , crosslinking 0j§], etc.) to perturb fluid 
phase separation, but the investigation of geometrical 
perturbations is more recent. There have been numer- 
ous studies on the perturbation of phase separation aris- 
ing from the presence of a plane wall, which is one of 
the simplest examples of a geometrical perturbation of 
phase separation |9[-|ll|. Measurements and simulations 
both show the development of "surface-directed" compo- 
sition waves away from plane boundaries under the condi- 
tion where one component has an affinity for the surface. 
The scale of these coarsening surface waves grows much 
like those of bulk phase separation patterns jBllq] • Re- 
cent simulations pj,[lq] and experiments Jf^M p2[ have 
shown that variation of the polymer-surface interaction 
within the plane of the film allows for the control of the 
local polymer composition in blends phase separating on 
these patterned substrates ( "pattern-directed phase sep- 
aration" |17|l ). Measurements have also indicated that 
the polymer-air boundary of phase-separating blend films 
on patterned substrates can be strongly perturbed by 
phase separation within the film p2|-p4f and thermal fluc- 
tuations of the polymer-air boundary can also strongly 
influence the structure of thin polymer films |25| . 

In the present paper we focus on the consequences of 
having geometrical heterogeneities of finite extent in a 
phase separating blend. The Cahn-Hilliard-Cook (CHC) 
theory [g(| for phase separation is adapted to describe 
phase separation of a blend with spherical, cylindrical 
(fiber), and plate-like shaped filler particles. The ex- 
tended dimensions of the fiber and platelet filler parti- 
cles are taken to be much larger than the scale of the 
phase separation process. A variable polymer-surface in- 
teraction is incorporated into the filler model in a fashion 
similar to previous treatments of plane surfaces [O [14J . 

The paper is organized as follows. In section II we 



briefly summarize the CHC model to introduce nota- 
tion, to define the relation between model parameters and 
those of polymer blends theory and to explain modifica- 
tions required for incorporating immobile filler particles 
into the CHC simulations of phase separation. Section 
III summarizes the results of simulations for representa- 
tive situations. Key phenomena are identified: 1) Target 
composition patterns form in near critical composition 
blends. 2) Target patterns are a transient phenomenon. 
3) The scale of the target patterns depends on quench 
depth and molecular weight. 4) Qualitative changes in 
the filler-induced composition patterns occur when the 
surface interaction is neutral and when the blend com- 
position is off-critical. 5) Multiple filler particles induce 
composition waves exhibiting complex interference pat- 
terns. Section IV provides a simple analytic estimate of 
the scale of the target pattern based on the linearized 
CHC theory and these results are tested against sim- 
ulations for circular filler particles. In section V simu- 
lation results are compared to experiments on ultrathin 
polymer films having silica bead filler particles immo- 
bilized by the solid substrate on which the films were 
cast. Atomic force microscopy measurements on the filled 
blend films are compared to the analytical predictions. 
Simulations of phase separation in off-critical blend films 
are briefly compared to analogous experiments. Measure- 
ments on crosslinkcd blend films are also considered. The 
final section (VI) discusses generalizations of the present 
study to manipulate the structure of phase-separating 
blends. 



II. THE MODEL 



A. CHC Equation 



We present a brief discussion of the CHC model 1 26 2J-J 
to introduce notation and to explain the modifications 
required for incorporating filler particles. The modeling 
of the phase separation dynamics is based on gradient 
flow of a conserved order parameter </>(r, £), 
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(1) 



with 0(r, t) equal to the local volume fraction of one of 
the blend components. Incompressibility of the mixture 
is assumed so that the local volume fraction of the second 
component is 1 — <j). The mobility M is assumed to be 
spatially uniform and independent of concentration and 
the free energy functional F[(j)] has the general form 
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where f(4>) is the bulk Helmholtz free energy per lattice 
site, v is the volume per lattice site, k B is Boltzmann's 



constant, T is temperature, and k((J>) is a measure of 
the energy required to create a gradient in concentration. 
Higher order gradient terms are neglected. The chemi- 
cal potential is given by \x eq = df/d/fil^^, which ensures 
that </>(r) = (j) eq is the solution of (5F[(/)]/(50(r) = 0. Fi- 
nally, thermal fluctuations necessary to ensure a Boltz- 
mann distribution of 4>{r) in equilibrium are included via 
the Gaussian random variable £. The average of £ van- 
ishes, and £ obeys the relation 

(C(r, tX(r',t')) = -2Mk B TX7 2 6(r - r')S(t - t 1 ). (3) 

For temperatures near the critical temperature T c the 
free energy can be expanded in powers of the composition 
fluctuation ip( r ) = </K r ) — </>c, giving the Ginzburg-Landau 
(GL) functional 

F[1>(t)] = k B T J ^[I Kc (V^) 2 + \c^ + ±W 4 + ■ ■ •] 

(4) 

where k c = k(4> c ). Neglected terms are either higher 
order in ip or in ^/c ex \JT — T c , which is small near the 
critical point. Eq. (fil), in combination with (0), defines 
the well-known Model B |29|, 
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Eq. (ph is used to study the dynamics following a quench 
to the two-phase region T < T c , where c < 0. In that 
case, the CHC equation may be rescaled into the dimen- 
sionless form iBOf 
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by making the substitutions 

r^(|c|/« c ) 1/2 r 

t -> (Mk B Tc 2 /vK c ) t 

V^Mc|) 1/2 V 

Note that this amounts to rescaling space by V2£-, 
where £_ is the thermal correlation length in the two- 
phase region, and time by r = D co u/(2^_), with D co u 
the collective diffusion coefficient. Here the noise term 
f]{r,t) satisfies (r)(r,t)) = and 

(7?(r, t) V (r',t')) = -V 2 S(r - r')5(t - t'). (8) 

The only parameters left to specify the dynamics are the 
(conserved) average concentration ipo = (4>) — 4>c and the 
dimcnsionlcss noise strength parameter e, 
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Roughly speaking, the reciprocal of e is a measure of 
quench depth. The parameter e also arises in discussions 
of the width of the critical region, and the connection be- 
tween thermal noise strength and the Ginzburg criterion 
was first noted by Binder {HI . 



B. Polymer Blends 



C. Surface Energetics 



For polymer blends we take the Flory-Huggins (FH) 
form of the Helmholtz free energy per lattice site, 



In the presence of a surface we add a local surface 
interaction energy to be integrated over the boundary, 
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Here \ represents the monomer-monomer interaction en- 
ergy, Ni is the polymerization index of component i, and 
4> is the volume fraction of component A. For the coeffi- 
cient of the gradient term we use de Gennes' random- 
phase approximation (RPA) result (neglecting the en- 
thalpic contribution pTJ), 
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where a a and gb are monomer sizes of the A and B blend 
components, given in terms of the radius of gyration of 
the ith component R g> i by of = B? gi /Ni. 
FH theory exhibits a critical point at 

c = < 2 /(< 2 + < 2 ) 

Xc = X(4>c T c ) = [N A /2 + N B /2 ] 2 /{2N A N B ). (12) 

Consequently, the coefficients of the GL functional are 
defined as J32T] 



The coupling constant h in the leading term plays the 
role of a surface field which breaks the symmetry be- 
tween the two phases, i.e., attracts one of the compo- 
nents to the filler surface. The coupling constant g in the 
second term is neutral regarding the phases, and results 
from the modification of the interaction energy due to 
the missing neighbors near the surface |33| and chain con- 
nectivity |3j| . For studies of surface critical phenomena 
p5| , p6[ | and surface dynamics |37],[58[ one typically keeps 
only these terms as a minimal model of phase separation 
with boundaries. 

There has been considerable attention given to the sub- 
ject of the appropriate dynamical equations for the sur- 
face boundary conditions p7\ . We follow most authors 
and impose zero flux at the boundary, n • j,/, = 0, which 
gives 



h-V(V 2 ip + ip-ip 3 ) = 0. 



(17) 



For the second condition we impose that of local equilib- 
rium at the surface, namely, 
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The phase separation dynamics of polymer blends can 
then be described by the dimensionless CHC equation 
with e determined by molecular parameters, 



While more sophisticated treatments are available jl6|, 
they lead to dynamics which rapidly relax to equilibrium 
and satisfy the above condition. The details of how we 
implement (|17| ) and (|l8| ) in simulations with a curved 
interface are presented in Appendix A. 
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where N A = N and x = N B /N, f(x) = (1 + ^/x) 3 /8x, 
and 7 is the ratio of monomer size to lattice size, 
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7 simplifies to 7 — <ja,b /v 1 ' 3 for symmetric blends. Thus 
we see that deep quenches (x 3> Xc) and high molecular 
weight polymers effectively reduce the thermal fluctua- 
tions in the rescaled dynamical equations. 



D. Simulation Details 

The equation of motion is solved using a standard cen- 
tral finite difference scheme for the spatial derivatives, 
and a first-order Euler integration of the time step |39| . 
In all the simulations, the lattice spacing is taken between 
0.7 and 1.0 in dimensionless units (sufficiently smaller 
than any relevant physical length scales) and the time 
step is taken sufficiently small to avoid numerical instabil- 
ity. In the present paper, all simulations are performed in 
d — 2 on lattices up to size 128 2 , depending on the choice 
of mesh size. We note that the CHC equation is known 
to exhibit quantitatively similar pattern formation and 
coarsening kinetics in two and three dimensions. Fur- 
ther technical details about this type of simulation can 
be found in Ref. [El. 



III. ILLUSTRATIVE SIMULATIONS 

In Fig. [j] we show the influence of an isolated circu- 
lar filler particle on the development of the phase sepa- 
ration pattern of a blend film having a critical composi- 
tion. The thermal noise is small in these two-dimensional 
simulations (e = 10~ 5 ), corresponding to high molecular 
weight and/or a deep quench (low and high temperatures 
relative to the critical temperature for upper and lower 
critical solution type phase diagrams, respectively). At 
an early stage of the phase separation process the filler 
particle creates a spherical composition wave disturbance 
that propagates a few "rings" into the phase separating 
medium in which the filler particle is embedded. The 
target rings initially have the size of the maximally un- 
stable "spinodal wavelength" Ao obtained from the lin- 
earized theory |26|] . As the characteristic scale of the bulk 
phase separation pattern coarsens to the size of the filler 
particle, the outer rings of the "target" pattern become 
disconnected and increasingly become absorbed into the 
background spinodal pattern. The perturbing influence 
of the particle becomes weak at a late stage of the phase 
separation where the scale of the background phase sep- 
aration pattern exceeds the filler particle size. The finite 
extent of the filler thus limits the development of the 
composition waves to a transient regime. 
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FIG. 1. CHC simulation of the influence of filler particles 
on polymer blend phase separation in a critical composition 
blend. Calculations are performed in d = 2 with e — 10 -5 , 
g — 1.0 and h = 1.0. In the reduced units of CHC theory 
[cf. Eq. (@)], R = 5.6 (or in explicit units R ~ 3.96£"). 
The target phase separation patterns are well developed at 
early times, but fragment as the "background" spinodal phase 
separation pattern coarsens to a scale larger than the filler 
particle (central gray region in figure). 



The formation of target patterns can be described as 
composition waves propagating into the bulk, unstable 
region until they are overwhelmed by the developing 
background spinodal decomposition pattern. The rate 
of onset of spinodal decomposition is controlled by the 
strength of thermal fluctuations, hence the noise param- 
eter e plays a crucial role in determining the radial ex- 
tent of the composition waves. Fig. shows two sys- 
tems with equal surface interaction, but with varying 
noise strengths: (a) e = 10~ 3 and (b) e = 10~ 5 . We 
see that the spatial extent of the target pattern is larger 
for smaller e. Consequently, we expect deeply quenched 
and/or high molecular weight polymer blends to be fa- 
vorable systems for observing filler induced composition 
waves because of the relatively low thermal noise level 
typical of these systems, and the relatively high viscosity 
of these fluids which slows the dynamics and makes mea- 
surements of intermediate stage patterns possible (e.g. 
via atomic force microscopy as discussed in section V). 
In section [V, we use the linearized CHC equation to 
estimate the extent of the composition wave. 




FIG. 2. Composition waves resulting from different values 
of thermal noise: (a) e = 10~ 3 and time t — 21 in dimen- 
sionless units; (b) e = 10~° and time t = 33. The surface 
interaction parameters are h — 1.0 and g — for both, and 
the filler radius is R = 5.6. 

We observe that the composition waves disappear 
when the particle radius R becomes vanishingly small, 
and the persistent waves developing from planar surfaces 
are recovered for very large spherical particles. We then 
examine particles of sizes intermediate between these ex- 
treme limits and during the intermediate phase separa- 
tion period in which the target patterns are well devel- 
oped. Fig. g shows the angular-averaged composition 
profile ijj(z), with z = r — R the radial distance from the 
surface of the filler particle. We observe that the ampli- 
tude of the local composition fluctuations become more 
developed and more sharply defined with increasing filler 
size, R. In comparison to the planar surface (R — > oo), 
the composition wave profile for R = 10 is only slightly 
reduced in amplitude, whereas for R = 3 the amplitude 
is reduced to half that of the wall case. We also remark 
that the radial extent of the target pattern is similar for 
all particle sizes. 

We next examine the influence of the surface interac- 
tion on the formation of filler- induced target patterns. 
The impact of the symmetry breaking perturbation of 



the filler particle on the phase separation may be tuned 
through the surface interaction parameters g and h. We 
focus our attention on h since it has a predominant effect 
on the resulting pattern formation. 
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FIG. 3. Typical angular- averaged radial composition pro- 
files ip{z) for R = 3, 10 and oo in an intermediate stage of 
phase separation. Parameters used are h — 0.1, g — 0, 
e = 10 -2 , and the time is t — 25. Averages were taken over 
ten independent configurations. 




(a) h=2.0 
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FIG. 4. Influence of surface interaction on filler-induced 
pattern formation. The figures show intermediate stage phase 
separation (t — 50) with thermal noise e = 10 -4 and ra- 
dius R — 9.6. Filler-induced target pattern for a selective 
polymer interaction (a) h — 2.0 and (b) h — 0.1, and (c) a 
non-selective polymer interaction, h = 0. Note the tendency 
for the domains to align perpendicularly to the neutral filler 
interface. 

Fig. [| shows CHC simulations of blend phase separa- 
tion for the case where one component strongly prefers 
the filler, weakly prefers the filler, and has no preference 
for the filler (h = 0). The quench depth parameter and 
computation times are identical in these images. We see 
that the target patterns do not form in the case of filler 
particles with a ( "neutral" ) non-selective interaction, but 
rather there is a tendency for the patterns to align per- 
pendicularly to the interface. This type of compositional 
alignment, which we find is even more pronounced in the 
case of a planar surface, also occurs in block copolymer 
fluids |4t{ . The perturbing influence of the boundary in- 
teraction saturates with an increase of h, as can be seen 
by comparison of the h = 2.0 and ft = 0.1 cases. Be- 



low we demonstrate that the extent of the target pattern 
depends logarithmically on h. 

Target waves are a variety of spinodal pattern with a 
symmetry set by the shape of the filler particle boundary. 
The introduction of surface patterns on a solid substrate 
can similarly break the symmetry of the phase separation 
process and can be used to impart a particular "shape" 
to the spinodal pattern |L7 20 1. A previous CHC simu- 
lation by us considered this "patterned-directed" phase 
separation in near critical composition blend films |L7|. 

The blend composition also can have a large influence 
on the character of the filler-induced phase separation 
structures in blends, particularly when sufficiently far off 
critical to suppress the spinodal instability. In this case 
we find a layer of composition enrichment ("encapsula- 
tion layer") forms about the filler particle, but there are 
no target patterns Q|. This encapsulation layer grows 
in time, but appears to grow slower than t 1 ' 3 . A non- 
selective interaction (h — 0) leads to the absence of en- 
capsulation by the minority phase, and minority phase 
nucleation occurs largely unaffected by the presence of 
the filler. We thus find that the development of target 
patterns requires the conditions of ordinary spinodal pat- 
tern formation (i.e., "near" critical composition) and the 
existence of a heterogeneity to initiate the wave distur- 
bance. 




FIG. 5. Illustrative CHC simulation of blend phase separa- 
tion with many filler particles. Filler particles preferential to 
each phase are included (the two filler types have a dark and 
light filler core). Simulation parameters are t = 80, h — ±1.0, 
R = 2.0, and e — 10 -6 . Note the interference pattern between 
these composition waves. 

A representative example of the interference between 
filler-induced rings at a non-vanishing filler concentra- 
tion is shown in Fig. pi Filler particles can each have an 
affinity for the different blend components so that the en- 
riching phase ( "charge" ) can vary near the surface of the 
filler particle at the core of the target waves. (Sackmann 
has noted that composition enrichment patterns occur 



about membrane proteins in lipid mixtures comprising 
living cell membranes, and these patterns mediate pro- 
tein interactions, leading to attractive or repulsive inter- 
actions depending on the "charge" J42].) In the low noise 
limit it should be possible to obtain novel wave patterns 
like those found in reaction diffusion models with regu- 
larly spaced sources for wave propagation J&j] , but we do 
not pursue this here. We do mention that the use of filler 
particles responsive to external fields could allow the ma- 
nipulation of the large scale phase separation pattern if 
the external fields are used to align the filler particles. 



IV. ESTIMATION OF THE SPATIAL EXTENT 
OF THE TARGET PATTERNS 

In this section we derive an approximate analytical ex- 
pression for the spatial extent of the composition wave 
pattern. Our method is based on the observation above 
that the composition wave propagates until it is over- 
whelmed by the growth of the bulk spinodal decompo- 
sition background pattern. In the context of surface- 
directed spinodal decomposition, a qualitative explana- 
tion of this type of phenomena was proposed by Ball and 
Essery |12| , who argued that the early time dynamics can 
be adequately described by the linearized CHC equation, 
in which the composition wave and the bulk spinodal 
decomposition add linearly. Both processes continue in- 
dependently until a local non-linear threshold value of 
\ip\ ~ ipt is reached, which then relaxes toward the equi- 
librium value ip = ±1. 

We test this conjecture numerically in the case of the 
filler inclusions. First, we simulate the CHC equation in 
bulk, with no filler particle, and determine the time to 
at which the root-mean-square concentration (^ 2 (i)) 1 ' 2 
reaches the threshold value ipt = 0-15 ( our reason for this 
choice is given below). Next, we simulate the filled blend 
in the absence of noise solving for the pattern at time to, 
and then estimate the radius rg at which the composi- 
tion wave (envelope) exceeds ipt . Finally, we perform the 
simulation with both the filler particle and the thermal 
noise, and find the pattern to be well-characterized by 
the size r^ for a range of noise and surface interaction 
parameters. 

Based on these ideas, we next develop an analytic es- 
timate of the spatial extent of the pattern using the lin- 
earized CHC theory [£6|. At early times the order pa- 
rameter does not deviate significantly from zero, and one 
linearizes (0) to obtain 
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where fcg 



1 — 3-00 • We apply this equation first to 
the determination of the root-mean-square order param- 
eter ip T ms(t) = (^(i) 2 ) 1 / 2 , which characterizes the rate 
of growth of the concentration fluctuations at very early 



times following a quench to the two-phase region. Fourier 
transformation gives 



dti>{k, t) = k 2 {k 2 - k 2 )^j + e 1/2 r?(k, t). 



(20) 



The structure factor S(k,t) = (ip(\s.,t)ip(—\s.,t)) is then, 

S ^= { 2(*g-*») • (21) 

V'rms is found by integrating S(k, t) with respect to k, 
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Following frL^I , we observe that the integrand is sharply 
peaked about k = fco/v2 and we approximate it by a 
Gaussian. The integral is then readily evaluated to find 
(to leading order in 1/i) 
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This result, when tested against simulations of the fully 
nonlinear CHC, agrees well up to {ip 2 ^)) 1 / 2 w (0.15), 
thus motivating our choice for ip t indicated above. 

Finally, we equate (ip 2 (to)) to ip 2 to determine the time 
to at which the bulk phase separation process has reached 
the nonlinear threshold. The resulting transcendental 
equation for to can be approximately solved by observ- 
ing that ip 2 /e 3> 1 in low noise conditions, and that this 
ratio must be compensated primarily by the e fc °*°/ 2 fac- 
tor. Equating these two and then iteratively improving 
the estimate of to yields 
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Next we consider the linearized theory for the filled 
blend. We solve ( |l9| ) for the exterior of the filler par- 
ticle or fiber, and in the absence of thermal noise (the 
noise can simply be averaged out of the composition wave 
within the linearized theory). We consider idealized filler 
particles that are symmetric and finite in some of their 
coordinates and infinite (i.e., very large on the phase sep- 
aration pattern scale) in the remaining coordinates defin- 
ing the particle dimensions. With such symmetry, the 
composition wave depends only on the coordinate per- 
pendicular to the interface, which is a radial coordinate 
in d± dimensions. For example, a spherical particle in 
d = 3 corresponds to d±_ =3, a cylindrical fiber is pre- 
scribed by d± = 2, and a platelet filler reduces to the 
planar surface with d± = 1. We can treat the general d± 
case through the dj_-dimensional Laplacian, 



V 2 ip(r) 



d 2 ip 
dr 2 
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2 Or' 



(25) 



yielding a fourth order partial differential equation for 
( |l9| ) . The "source" for the composition wave comes from 
the boundary conditions obtained by linearization of (Jlj 
and (Jig), namely 



(26) 



f • V(fc2 + V 2 )iJj(R) = 

for conservation at the boundary, and 

f-Vip(R) = h + gip(R). 



(27) 



This provides a pseudo one-dimensional system which 
can be readily integrated numerically. 
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FIG. 6. Influence of particle shape on the development of 
filler-directed composition waves at early time. Radial com- 
position profiles ip(z) are obtained from the linearized CHC 
equation for symmetric particles finite in d± directions, cor- 
responding to "platelet" fillers (dj_ = 1) , fibers (dj_ = 2) , and 
spherical fillers (dj_ = 3). The radius is R — 10, and the 
surface interaction parameters are h = 0.005 and g — 0. The 
reduced time is t = 25. The inset shows the approximate solu- 
tions of (E9[) and (pll) (dashed lines) , compared to the d± = 1 
and 3 numerical solutions from the main figure (solid lines). 

The solutions for dj_ = 1,2,3 presented in Fig. g illus- 
trate the influence of d± on the composition wave pat- 
tern. Increasing d± reduces the amplitude of the com- 
position wave. This feature can be understood to arise 
from the increasing volume occupied by the outer rings. 
The opposite situation should hold for exterior bound- 
aries having these symmetries, so that more coherent ring 
structures might be anticipated in phase separation con- 
fined to these geometries, especially for spherical cavities. 
It may prove interesting to examine phase separation in 
the presence of fractal filler particles (like fumed silica) 
to determine whether geometry stabilizes or destabilizes 
the phase separation pattern and how the evolving phase 
separation pattern accommodates the fractal boundary 
structure. 

Returning to the analytical estimation of the target 
pattern size, we expand tp(r) in a basis which diagonal- 



izes the Laplacian. This amounts to performing a co- 
sine transform for d± = 1, a (Jo) Hankel transform in 
d = 2, and a half-integer Hankel tranform in d = 3. The 
last can be re-expressed as a Fourier cosine transform of 
rtp(r) rather than ip. Here we study the extremal cases 
of di = 1 and 3, and show that d± has only a minor 
effect on the pattern size. 

First, we revisit the d± = 1 case already addressed by 
Ball and Essery [Eg . One can solve (03) with boundary 
conditions via Fourier cosine transformation with respect 
to r and Laplace transform with respect to t, with the 
result, 



4>(k,s) 



hk 2 



s[s-k 2 {kl-k 2 )] 



(28) 



for the case g = (to which we restrict our attention). 
Inverting the Laplace transform and using a Gaussian ap- 
proximation again to invert the Fourier cosine transform 
gives the solution 



Vw=i(M)«V>(0,i)e- (z/4Ao)2 cos(z/A ), 



(29) 



^ v /^exp(^), Aq- 1 - fc /V2, and 



where V(0,i) ee ^Aexp^iJ, A" 1 

where subdominant terms in 1/t have been neglected. 
A non-zero value of g, while complicating (Eq), would 
appear in this approximate solution only via the sub- 
stitution h —* h + gtpo- hi this d± — 1 example, the 
R dependence drops out of the linearized equation, and 
z = r measures the distance from the wall. 

For a given noise strength e we have a time to at which 
the local composition of one phase reaches ip t . For sur- 
face interaction h we solve for the distance z out to 
which the envelope of ip(r, to) exceeds ipt. This gives the 
approximation, 



Zq 



2k Q tQ 




(30) 



Thus, the propagation front grows with a velocity 2fcg at 
long times [L2J , and the terms in square brackets are the 
leading correction to this long time asymptotic behavior. 
The d± = 3 composition profile may be obtained by 
observing that ( |19|) and ( |25|) yield the same equation for 
ripd ± =3(z, t) as for the composition profile ipdj_=i(z,t). 
Thus, we impose the boundary conditions (to leading or- 
der in R/r) and follow the above derivation with the 
result 



ipdj_=3(z,t) 



R 



R + z 



ip d±=1 (z,t). 



(31) 



Solving for the value zq at which the composition wave 
envelope equals ipt leads to nearly the same expression 
as the d± = 1 case (|30|), with an additional — ln(l + 
Zo/R)/kobo term in the square brackets. Typically zq < 



lOi?, fco is of order unity, and to ranges from 10 to 30, 
making this term roughly a 10% correction. 

To compare with our simulations, we consider d± = 
d = 2. For zq we simply take the arithmetic mean of the 
values obtained for d± = 1 and d± = 3 (motivated by 
numerical solutions). Equation (|24|) for t is substituted 
into (pOI) to obtain a prediction for Zq in terms of h, e, 
and 7/V Here we consider critical quenches with ipQ = 0, 
or ko = 1. Finally, we approximate lnio and lnzo with 
typical values, which introduces less than 10% error with 
the range of parameters considered here, and thus obtain 



za 



-2.5 ~ 8.5 log 10 e + 4.6 log 10 h. 



(32) 



A similar expression results for d = d± =3, with the 
primary difference being a change of the log 10 e coefficient 
to -7.9. 
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FIG. 7. Comparison between analytic estimate of the spa- 
tial extent of the target pattern the corresponding CHC sim- 
ulation at an intermediate stage of phase separation. Sim- 
ulations are performed at noise levels (a) e = 10~ 2 corre- 
sponding to modest molecular weight and shallow tempera- 
ture quenches, (b) e = 10 -4 , and (c) e = 10 -6 , correspond- 
ing to high molecular weight mixtures and deep temperature 
quenches. The surface interaction h — 1 for all cases, the 
radius R = 5.6, and the times for each quench are chosen to 
correspond to when the spinodal pattern is fully developed: 
(a) t = 20, (b) t = 30, and (c) t = 40. The solid lines are the 
predictions of (p3). 

Fig. shows a comparison of this estimate with the 
simulations. We see that the analytic approximation pro- 
vides a good rough estimate of the spatial extent of the 
phase separation pattern, although it predicts a size typ- 
ically one oscillation larger than the outermost unbroken 
target. 

We point out that the spherical composition waves are 
apparent in the average composition profiles even in the 
rather noisy looking ring patterns found in the late stage 
of target pattern formation. In Fig. pi we show a target 
pattern at intermediate values of noise, as well as the ra- 
dial average of the composition profile about the center 
of the target. Comparison shows that the ring composi- 
tion pattern persists in the radial average even after the 
target pattern appears visually to have broken up. This 
provides a possible explanation of the apparent overesti- 
mate of the target size in Fig. 7. 




(b) 




FIG. 8. Late stage target pattern and corresponding ra- 
dial composition profile, a) Filler-induced phase separation 
pattern for an intermediate noise value of e = 10 -3 at time 
t = 35, h = 1, and R — 5.6. b) Averaged radial profile about 
center of filler particle. Note that the ring composition pat- 
tern persists in the radial average even after the target pattern 
appears visually to have broken up. 



V. COMPARISON WITH EXPERIMENTAL 
RESULTS 

While scattering measurements of the growth of com- 
position waves are readily performed for a single plane 
boundary |l0[ , these measurements become more difficult 
in filled blends where the filler particles are randomly dis- 
tributed within the blend. This situation is unfortunate, 
given the predicted transient nature of the composition 
wave patterns when the particles are small. However, 
real space studies of blend phase separation are possi- 
ble in films sufficiently thin to suppress the formation 
of surface-directed waves normal to the solid substrate 
p4| . Under the favorable circumstances that one of the 
polymer components segregates to both the solid sub- 
strate and the polymer-air boundary of these nearly two- 
dimensional ( "ultrathin" ) blend films, phase separation 
is observed within the plane of the film |23|,g4],[44|] . The 
variation of the surface tension in the film accompanying 
phase separation gives rise to film boundary undulations 
that can be measured by atomic force microscopy (AFM) 
and optical microscopy (OM) |2J,Q. The thickness of 
ultrathin blend films is typically restricted to small val- 
ues (L rs 200 nm) and the height contrast of the surface 
patterns tends to become larger in still thinner films p2[ . 
A film thickness in the range of 20 - 50 nm is often suited 
for observing well resolved phase separation surface pat- 
terns similar to those found in simulations of bulk blends. 
In the following we compare our results with those of a 
model blend utilized in ultrathin phase separation studies 
reported elsewhere with silica beads added as the model 
filler J45|. 

The spun-cast films are composed of a near-critical 
composition blend of polystyrene and poly(vinyl methyl) 
ether (PVME). The filler particles are silica beads having 
an average size of about 100 nm, as measured by direct 
imaging of the particles. This particular filler was chosen 
because of its tendency to be enriched by polystyrene, 
rather than PVME, which enriches both the solid and 



air surfaces. In this way, the filler particles are not com- 
peting with the solid or air surfaces for the enriching 
polymer. Phase separation was achieved by annealing 
the film approximately 15° within the two-phase region, 
corresponding to a fairly shallow quench. Film topogra- 
phy (height) was measure by AFM. Further details of the 
experiment are provided in Rcf. p5[. 




in the radially averaged patterns and in Fig. |l(J we show 
the radial average of the AFM height data centered about 
a representative filler particle. The target pattern in the 
radially averaged data extends far beyond the ring fea- 
ture apparent in the image in Fig. EH The data in Fig. HG 
correspond to a shallow quench, and are comparable to 
the intermediate stage, shallow quench simulation data 
in Fig. | 

Next we directly compare the prediction of the lin- 
earized theory to the AFM data. An exact solution of 
ipd ± =2(z 1 t) is difficult, but we can obtain a reasonable 
approximation to ipd ± =2{z, t) by generalizing the method 
described above for tpdj_=i(z, t). We estimate ipd±=2{z, t) 
as a Gaussian decay function multiplied by the eigenfunc- 
tion of the Laplacian in d = 2 [rather than in d = 1 as in 
the case of (|29[)]. In this approximation, ipd ± =2(z,t) be- 
comes a product of a Gaussian as in Eq. (|29| ) and a Bessel 
function J (2irz/\ ), and we show a fit of this function to 
the AFM data in Fig. [l(]. The fitted value of the particle 
radius R is 82 nm, which is comparable to the average 
particle radius obtained by optical microscopy (R sw 100 
nm). The scale parameters of the Gaussian and Bessel 
functions have been adjusted along with the prefactor 
which is set by the value of i/jdj_ = 2(z,t) as z tends to 
zero. It is clear that the oscillatory pattern scale is on 
the order of the background phase separation pattern, 
and that the linearized expression for ipd ± =2{z> t) has the 
qualitative shape of the measured profile. Such quali- 
tative agreement is the best that can be expected from 
the linearized theory, which strictly speaking should hold 
only at very early times. 



FIG. 9. AFM image of phase separation pattern in 
PS/PVME blend film with dilute concentration of filler parti- 
cles at (top) 20/irn and (bottom) 100 /im scale. The height un- 
dulations reflect composition variations within the film asso- 
ciated with surface tension variations J44J . Image contrast has 
been enhanced by a film washing procedure. From Ref. []45|. 

Fig. |9|(top) shows the topography of the blend film at 
an intermediate stage of phase separation where we ex- 
pect circular filler-induced composition waves to be evi- 
dent. The pattern resembles the simulated patterns un- 
der similar quench conditions. The symmetry of the film 
phase separation pattern is locally broken by the presence 
of the filler particles, leading to the formation of ring-like 
concentration wave patterns. Note that when observed 
on a larger scale [Fig. ^(bottom)], the phase separation 
pattern far from any filler particles resembles the typi- 
cal spinodal decomposition pattern observed in control 
measurements on the same blends without filler. 

It is apparent that the patterns in Fig. are in a rel- 
atively late stage of phase separation, where the rings 
are beginning to break up along with the "background" 
phase separation pattern. The simulations above indicate 
that the target patterns are more persistently expressed 
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FIG. 10. Radial average of a two-dimensional fast Fourier 
transform of AFM height data centered about an isolated tar- 
get pattern in Fig. H. The scale of the damped oscillatory 
profile is reduced by the phase separation scale determined 
by a Fourier analysis of the AFM height data for the entire 
film. The solid line is a fit to the data as described above. 
Data from Ref. IBI. 



At still longer times, the phase separation pattern 
eventually breaks up into droplets and little difference is 
observed between the films with and without filler. Thus, 
the target patterns induced by the filler particles are tran- 
sient, as observed in the simulations. Of course, the ver- 
sion of the CHC model used here cannot reliably describe 
quantitative features of these late stage processes without 
the incorporation of hydrodynamic interactions. 

Under far off-critical conditions and a selective interac- 
tion between the filler particles and one of the polymers 
(h > 0) the filler particles are "encapsulated" by a layer 
of the favored polymer so that concentration waves do 
not develop. The formation of droplets by nucleation or 
far off-critical spinodal decomposition can also have the 
effect of breaking the symmetry of the phase separation 
process, but the pattern formation is not generally the 
same as for critical composition mixtures. Recent mea- 
surements have reported the occurrence of filler encapsu- 
lation in a blend of polypropylene and polyamine-6 with 
glass bead filler particles |§,0- Encapsulation occurs 
when the polypropylene-rich phase having the selective 
interaction for the filler is the minority phase, but no en- 
capsulation occurs when polypropylene is the majority 
phase. This finding compares well with the simulation 
results discussed in section IV. 




FIG. 11. Phase contrast microscopy image of phase sepa- 
ration pattern of a photo-crosslinked blend. The scale of the 
image is 10 /im. From [p|. Reproduced with permission from 
Marcel Dekker. 

Radiation crosslinking provides another source of het- 
erogeneity that can be introduced readily in phase sep- 
arating films. Measurements of irradiated photoreac- 
tive blends of PVME and PS with a crosslinkable side 
group styrene-chloromethyl styrene random copolymer 
(PSCMS) show the formation of striking ring composi- 
tion patterns p6|,(l7f and we reproduce one of these pat- 



terns in Fig. [n] (compare with Fig. ||). Furukawa J4q ] 
has interpreted these observations in terms of a model 
by which irradiation first brings the blend into the nucle- 
ation regime where droplets phase separate, followed by 
the entrance into the spinodal regime where the droplets 
act like the filler particles discussed in the present paper. 
This is a plausible interpretation of the qualitative origin 
of these patterns, but it is difficult to interpret these mea- 
surements directly from CHC simulations since crosslink- 
ing imparts a non-trivial viscoelasticity to the polymer 
blend EjJ. The crosslinking, which also increases the 
molecular weight of PSCMS, and the increased elastic- 
ity, both lead us to expect a decrease in the thermal 
noise and thus an increased tendency to form target pat- 
terns. It also seems plausible that the crosslinks them- 
selves provide the source of heterogeneity, inducing the 
development of composition waves. 

VI. DISCUSSION 

The presence of filler particles in a phase separating 
fluid mixture is found to give rise to transient composi- 
tion wave patterns in simulations based on the Cahn- 
Hilliard-Cook model and in measurements J45{ on ul- 
trathin polystyrene/poly(vinylmethyl) ether blend films 
with silica filler particles. In both the simulation and the 
experiment, the composition wave patterns were found 
to be transient and the filler is found to have a dimin- 
ishing effect as the scale of the phase separation pattern 
becomes larger than the filler particles. The propaga- 
tion of composition waves is enhanced at lower thermal 
noise level so that the effect propagates to larger dis- 
tances for deeper quenches and higher molecular weight 
blends. The finite size and the dimensionality of the filler 
particles are found in our simulations to have a similar ef- 
fect in determining the stability of the composition wave 
pattern at intermediate times. The composition waves 
become more stable for particles large in comparison to 
the spinodal wavelength and the concentration waves ex- 
hibited by these larger particles arc similar to planar in- 
terfaces. The composition waves about the filler particles 
are more stable for particles extended at great distance 
along more directions; i.e., surfaces are more stable than 
long cylinders, which are more stable than spherical filler 
particles. Our results compare favorably with experi- 
ments on phase-separating filled blend films which are 
nearly two-dimensional. 

Filler particles are an example of a perturbation of 
phase separation by boundaries interior to the fluid. It 
would be interesting to investigate the influence of ex- 
terior boundaries of finite extent on phase separation. 
It seems likely that composition waves within confined 
geometries should be more stable because of the decreas- 
ing surface area of the rings farther from the surface. 
This should lead to well developed and more long last- 
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ing perturbations of the phase separation process. The 
relation between boundary shape and phase separation 
morphology should be very interesting for this class of 
measurements. Phase separation within arrays of filler 
particles where the distinction between interior and ex- 
terior boundaries becomes blurred and where larger per- 
turbations of the phase separation process may be antic- 
ipated, should also prove interesting. The distinction be- 
tween large and small and fixed and mobile filler particles 
should lead to a range of new phase separation morpholo- 
gies since the development of composition waves should 
lead to changes in the filler-filler interaction that can in- 
fluence the subsequent development of the film structure. 
The utilization of geometrically and chemically patterned 
surfaces and additives offer many opportunities for the 
control of the phase separation morphology and resulting 
properties of blend films, and the study of these surface- 
induced phase separation processes raise many interest- 
ing problems of fundamental and practical interest. 
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APPENDIX A: BOUNDARY CONDITIONS ON A 
CURVED SURFACE 

Curved boundaries complicate the implementation of 
boundary conditions in a spatially discretized simulation. 
In the present work, we use a square lattice and simu- 
late filler particles with circular, cylindrical, and spher- 
ical symmetry. This requires a method of incorporating 
the boundary conditions that minimizes the effects of er- 
rors caused by approximating curved boundaries by lat- 
tices. In this appendix we present our approach to this 
problem. 

Generally, ( p7| ) and (18) are imposed by inclusion of ip 
and \jl = — V 2 ip — -0 + Rvalues at the lattice sites on the 
immediate interior of the boundary (within the wall or 
filler), which are determined from the boundary condi- 
tions before each time step. We superimpose the circular 



In both cases of Fig. [L2] we use the three vertices 
shown as open circles for the discretized representation 
of ip(xo,yo) and its normal derivative f ■ V^>(xo,yo)- 
To highest order these representations are unique [to 
0(Ax 2 ) for ip and O(Ax) for the derivative, with Ax the 
lattice spacing]. Hence, ([18]) may be used to determine 
ipi,j, the field at the interior site, from the appropriate 
exterior points. We find for case (a) the relation 



4>i. 



(1 — gtj [(sin 9 — cos 9) ipi.j+i + cos 9 ip. 



i+ij+ij 



hAx 



(1 - g£)sin9 + gAx 



(Al) 



where I is the distance from the interior lattice site to 
(xo,yo) while 9 is the angle between the radius and the 
horizontal axis. For case (b) the analogous expression is 



V 1 *,; 



(1 — g£) \cos6 V'i+iJ + shr^ V'ij'+i] — hAx 
(1 - .9^) (cos 9 + sin 9) + gAx 



(A2) 



The chemical potential fiij may also be assigned at the 
interior point, in practice by assigning a value to (V 2 ?/>)ij 
to supplement the Laplacian derived from ip outside the 
boundary. In this way we can impose the conservation 
requirement ([17]) for case (a) via 

Hij = (1 - cot 6) fM,j+i + cot 9 /i i+ i j+i , (A3) 

while for case (b), 



Mi,j 



■ Mi+lj + 



sin# 



IM,j+i- (A4) 



In simulations with thermal noise we assume a separation 
of time scales between thermal fluctuations and order 
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parameter variations (as described in |14|), and simply 
supplement the above conditions with the conservation 
law for fluctuations at the boundary: f • v = 0, where v 
is the noise current derived from 77 = V ■ v. 
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